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1. INTRODUCTION 



Finite state space semi Markov processes find application in many areas 

including queueing, reliability, and clinical trials, (cf Cinlar [1975], Weiss 

§ 

and Zelen [1965], Cox [1985]). Often interest in the application centers on 
the distribution of a first passage time to a state or a set of states 
representing for example the lifetime of a system or the end of a busy period 
of a server . 

In this paper we consider the problem of the estimation of the log 

survivor function of a first passage time for a fixed time t. We suppose that 
observational data are known about the semi Markov process in question. 

Although the approach and results are given concretely for one particular 
semi -Markov process, they apply more widely. 

To be specific, consider N individuals. Let X^(i) be the state of ith 
individual at time t. We will assume {X^(i);t>0} are independent identically 

distributed processes having the same probability law as { X ^ ; t > 0 > . The 
process {X^;t>0} is a semi-Markov process with three states {0,1,2}. State 0 
is absorbing and X^=l . The sojourn time in state i has a distribution 

function F . . Upon leaving state 1, the process transitions to state 0 with 

probability 9 and to state 2 with probability 1 - 9 . From state 2 the process 
transitions to state 1 with probability 1. 

Let 

D inf { t > 0 : X t - 0} 

the entrance time to state 0. 
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The problem is to estimate? the? logarithm of P{D>t}, for fixed t, from data 
obtained by observing the N individuals. Three possible estimation procedures 
will be considered. The 1 three estimators use different amounts of information 
concerning the process. One procedure uses only the observed entrance time to 
state 0 for the N individuals to estimate the empirical distribution function 
of D. Another procedure makes parametric assumptions concerning the sojourn 
time distributions and uses maximum likelihood to estimate the probability. 
A third approach uses an exponential approximation to P{D>t} and empirical 
distributions to estimate . 

In Section 2, the estimation procedure based on an exponential 
approximation to P{D>t} is described and asymptotic results are obtained for 
it . 

In Section 3 the other two estimation procedures are described and the 
results of a simulation experiment are given. The experiment was conducted to 
compare the performance of confidence intervals for the three estimators for 
small and moderate numbers of individuals. It is found that the procedure 
based on the maximum likelihood is the most efficient if the underlying model 
is correct. However, maximum likelihood confidence interval coverage is 
sensitive to incorrect model assumptions. The nonparamet r ic confidence 
intervals based on the exponential approximation have the correct coverage as 
long as the time t is not too small. The conservative binomial confidence 
intervals based only on the observed times of absorption are the least 
efficient and tend to overcover. 

2. AN ASYMPTOTIC RENEWAL ESTIMATE 

The probability P{D>t} for the model described in the Introduction 
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satisfies the equation 



t 

P{D>t} = g(t) + (1-0) S Q (F *F )(ds) P {D> t s } (2.1) 



with 



t 

g(t) - F 1 ( t ) + (IS) / o F^ds) F 2 (t-s). (2.2) 



This is a renewal equation with defective inter-renewal distribution 



L(t) - (l-^)(F l *F 2 )(t). (2.3) 

Following Feller [1966] page 326, assume there exists k which satisfies 
the following equation 

S°° e Kt L ( d t ) = (1-0) ♦ (k) * (k) = 1 (2.4) 

0 J “ 



where 



♦ .(C ) = e ct F 1 ( d t ) 



(2.5) 



If (F *F ) is not arithmetic, then under certain integrab i 1 i ty conditions 

lim e Kt P{ D>t ) = - (2.6) 

t-eo 
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whore 



OO V C 

fJ ~ / se L(ds) 
0 



(2.7) 



and 



b = /°° e Kb g(s)ds=^(K). 

_ K 1 



( 2 . 8 ) 



Following an approach that was found useful in another context, (Gaver 

and Jacobs ([1986]), an estimator for P{D>t} can be obtained by using 

estimates of k , b and fJ in (2.6). Specifically, let {S (i),n=l,..., M.} be 

n i 

the collection of sojourn times in state i for the N individuals and let R be 
the number of transitions from state 1 to state 2 for all the N individuals. 



Put 9 



N 

N*-R 



and 






M. 

l 



r 1 e tS k (i) 
k=l 



i = 1,2. 



(2.9) 



An estimator of k is the solution, k, to the equation 



( 1 -*) ♦ 1 (C)^ 2 («) = 1. 



( 2 . 10 ) 



Notice that the left hand side of equation (3.10) is nondecreasing in £ and at 
£ = 0 equals (1-9) < 1. Hence (2.10) has a unique solution which can be found 
numerical ly . 
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The estimators for the other quantities in (2.6) are 



„ . M 1 M 2 

H = (1 6) mV _ Z z ex P{« < S k < 1 ) +S J ( 2 ) )> (S k (l) * S j ( 2 ) ) (2.11) 



1 2 



k= l j--l 



b = t ^(k) . 

K 



( 2 . 12 ) 



The asymptotic exponential estimator for P{D>t} is 



P A (D>t) = £e' Kt . 



(2.13) 



2 . 1 Large Sample Properties of k and In P ^ (D>t) 



Let 

f(f.a) 1 (l-a)^ 1 ({)^ 2 ({) = 1 -(l-a)h(e). 

Expand h(x) in Taylor series about the solution k of (2.4). Since f(x,0) = 0, 
0 = f(K.e) = 1 (1-0) [ h (« ) + (k -k ) h ' (k ) + I (k K) 2 h"(a/t)] (2.14) 

= f(K,0) + (0 -0 ) h (/< ) ( 1-0 ) (K-K ) [h* (k ) + i(x -K)h* *(a*)]. 



Thus 



K K 



f (k ,0 ) + (0-0 ) h(K ) 



B 



(2.15) 
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where 



B = (1 0)[h‘(K) + h ,, (<ut)] 1 



(2.16) 



and a is a constant and 9 is the probability of a transition from state 1 to 
state 0. By the strong law of large numbers, as N-*» f(K,0)^O and 9 9-* 0. 
Further, as N-*», with probability one 

h* (*) -* (1-®) r se " S (F *F )(ds) (2.17) 

0 1 * 

and 

h(K) -► ^ («) i> 2 (K) . (2.18) 

Therefore, k-k-+ 0 as N-*», with probability 1 so the estimator k is consistent. 

To obtain a central limit theorem for k, the following moment condition 
is needed. Assume 



r e 2 * 55 F *F (ds) < «. 
0 



(2.19) 



Since f ( k ,$ ) = 0 , 



yT ( f(x.fl)] = /U [ f(«.e) f(*,«)] 



N (1-9) l* 1 («)* 2 («) ♦ 1 («)* 2 (k)1 



N (1-0) {[♦j(k) - ♦ 2 (K ) ] * 2 (/c) + * 1 (K)l* 2 (K)-t> 2 {K)]) ■ (2.20) 
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Thus, (2.19) and the central limit theorem imply that /If (f(K.0)] is 
asymptotically normal as N-*» with mean 0 and variance 

x T k T 

Var f = (1 9) e { ( 1 e )# 2 (K ) 2 Varfe ] ] + ♦ (k) 2 Var(e 2 ]} (2.21) 

where T. is a random variable with distribution F . 
i i 

Further the central limit theorem implies that as N -■», J ~ N [0-0] is 

2 2 3 

asymptotically normal with mean 0 and variance V. = 9 (1-0) * 0 (1-0). 

b 

Hence, as N-*°, (x x ) is asymptot i ca 1 1 y normal with mean 0 and variance 

V k " (V f * * 2 (*)) 2 V b )[(l-9) E[(T 1 -T 2 )e A:(T r T 2 ) ]]" 2 (2.22) 

A similar argument shows that i n P^{D>t} is also asymptotically normally 
distributed as N-^. 

3. S IMl'LATI ON RESULTS 

A simulation experiment was carried out to investigate small and moderate 
sample size behavior of confidence interval procedures for three estimators 
of £n P{D>t}. The following three estimators use different amounts of 

information about the semi-Markov processes { X ( i ) ; t>0}, i = l N, for the 

N individuals 

3.1 The Binomial Estimator 

An estimator which uses only the observed entrance times into state 0 for 

the N individuals, d , c! , . . ,d k ,. is 

1 2 N 
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P R (D>t> 



_ 1 _ 

N 



N 



(3.1) 



1 

1=1 



(t,< 



>)<V 



where 



if d . > t , 

l 

otherwise . 



The distribution of the estimator P {D>t} is binomial and a binomial 

B 

confidence interval [L ,IJ J can be constructed for P { D> t } cf. Larson [1982] p 

B B 

397. The binomial confidence interval for <n P{D>t} is (#n L_,fn U ). 

B B 



1 ( t ,00) ( d i 



3.2 The Asymptotic Exponential Estimator. 

The asymptotic exponential estimator for P{D>t} is given by P A (D>t} of 

(2.13). The estimator for the £n P{D>t) is in P.{D>t) and is asymptotically 

A 

normal as N-*». However, since the simulation experiment is for small to 
moderate numbers of individuals, two nonparamet r i c methods are used to obtain 
conf idence j ntervals . 

The jackknife is a procedure originally introduced by Quenouille [1956] 
for bias reduction, and adapted by Tukey [1958] to obtain approximate 

confidence intervals. Suppose interest is on a parameter p that is estimated 

by p using a complex calculation from data. The idea is that of assessing 

variability by recomputing p after removing independent subgroups of data and 
then using the recomputed p values to estimate a variance which is in turn 
applied to state a two-sided confidence interval that contains the true p with 
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specified confidence; see Efron [1982J and his more recent work, or Mosteller 
and Tukey [1977] for more details. For simulations involving N-20 
individuals, the jth subgroup consists of all data corresponding to the jth 
individual. For simulations involving N^50 individuals the first subgroup 
consists of all data corresponding to the first five individuals the second 
subgroup of all data corresponding to the second five individuals, etc. Some 
simulations for N = 50 individuals were run with the jth subgroup consisting of 
the data for the jth individual alone; according to some work these should be 
an improvement on those with 5 in a subgroup. However, the resulting 
confidence intervals differed little from those obtained by leaving out 5 
individuals at a time. 

The bootstrap is an alternate method that may be used to assess uncertainty 
of estimators; cf. Efron [1982] and his more recent work. In the simulations 

reported here a bootstrap replication is generated as follows. Let be the 
number of transitions from state 1 to state 2 for individual n. Let {S (i)} 

K 

be the collection of all sojourn times in state i for all individuals. To 

generate bootstrap data for one individual, one observation is drawn at random 

with replacement from {R }; call the observation r; r+1 observations are drawn 

n 

at random with replacement from the collection of state-1 sojourn times 

{S^(l)>; r observations are drawn at random with replacement from the 

collection of state 2 sojourn times {S (2)}. Bootstrap data is generated for 

k 

N individuals and the estimate £n P A (D>t) is computed using this data. This 
completes one bootstrap replication. B = 100 bootstrap replications are done. 
The B = 100 bootstrap estimates of ^n P A (D>t} are ordered. A 100 x (1 <t) % 
confidence interval is constructed using the cl/2 and l-a/2 quantiles of the 
bootstrap estimates. 
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3.3 The Max imum_ Likelihood Estimator 



In what follows, we assume the sojourn time in state i has an exponential 

1 



distribution with mean 






The log likelihood function under these 



assumptions is 



N 

L = l [R *n (1-0 ) ♦ €n 0 ♦ R *n p ♦ ( 1 +R )€n p (3.2) 

n n 2 n 1 

n=l 

V„ ( " W 2 ' 1 

where R is the number of transitions from state 1 to state 2 for individual n 
ri 

and T^(i) is the total time spent in state i by individual n. The maximum 
likelihood estimators are 



N+-R 
~ _ N+R 

p i " fTT) 

* = R 

P 2 " Tl2) 

N N 

where R = Z R , T(l) = Z T (1), and T(2) 
n n 

n = 1 n-1 



N 

Z T (2) . 
n 

n = l 



(3.3) 

(3.4) 

(3.5) 



Let D be the entrance time to state 0 for an individual; then 



V P 2 V p 2 6P 1 

S(t) = P(D>t } = ( n — - exp (A t ) - -L—t exp{A t } } j- (3.6) 

" 2 2 A 1 1 A l" A 2 
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where are t ^ le root:s °f the equation 

«PjP 2 f y(Pj + P 2 ) + y 2 = o. ( a . 7 ) 

The maximum likelihood estimator for P{D>t} can be obtained from (3.6) and 
(3.7) by replacing p^ , p^ and 9 by their maximum likelihood estimators and 
taking the logarithm of the resulting estimator of p{D>t). The distribution 
of the resulting estimator is asymptotically normal as tU°o. Since the maximum 
likelihood estimators 9 ,p ^ , and p^ are uncorrelated, the asymptotic variance 
i s 



Var[tn P{D>t}] = 1 [Var(e) (||-) + Var(p )(^L) 4 Var (p ) (£^_) ]. (3.8) 

09 l ap 2 op 

set r 

Confidence intervals are constructed using the asymptotic normality of the 
maximum likelihood estimator and the asymptotic variance with parameter values 
being replaced by their corresponding estimates. Confidence intervals based 
on the maximum likelihood ratio test may also be constructed [cf. Cox and 
Hinkley (1974) pp. 343] While these intervals are asymptotically equivalent 

to the normal confidence intervals that are computed, they may be better for 
sample sizes. 

3 . 4 The Simulation Experiment 



All simulations are carried out on an IBM 3033AP computer at the Naval 
Postgraduate School using the LLRANDOMII random number generating package; [cf 



Lewis and Uribe (1981)]. Some details of the simulation are given below. A 
more complete account can be found in Kim [1987]. In each replication 

estimates and confidence intervals are computed for the logarithm of P { I)> t } 
using the procedures described in subsections 3.1 -3.3. Each simulation is 
replicated 300 times. For each procedure, the number of confidence intervals 
covering the true value of £n P { D> t } is recorded. Also recorded are the 
number of intervals that arc too low (true ^ n P { D> t } > U, the upper endpoint 
of the interval) and too high (true £n P{D>t) < L, the lower endpoint of the 
interval). The average length of the confidence intervals is also computed. 

Tables 18 present results from simulations for N = 20 and N^50 
individuals. The first row in each table lists the confidence interval 
procedures: the binomial confidence interval for logarithm of estimate (3.1) 

(BIN); the maximum likelihood interval (MLE); the jackknife confidence 
interval for the asymptotic estimate, (JK); and the bootstrap confidence 
interval for the asymptotic estimate, (BT). Nomina] 80% and 90% confidence 
interval for £ n P{D>t) are computed using each procedure for each time 
t = 0.5, 1.0, 1.5, 2, 3, 4. Confidence limits that are greater than 0 are set 
equal to 0. 

Tables 1-4 present results for the case in which the sojourn times in 
state 1 have an exponential distribution with p ^ = 1 and the sojourn times in 
state 2 have an exponential distribution with mean p 1 = 0.1. The values of 
P{D>t} for the model for t = 0.5, 1.0, 1.5, 2, 3, 4 are given in Appendix A. 

Tables 1-2 show results for simulations with N=20 individuals. 

Table 1 shows coverage results. Under each procedure is given the number of 
intervals that cover the true value (C), the number of intervals that are too 
low (L) and the number of intervals that are too high (H) . In parentheses 
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next to these numbers are given the corresponding fractions. If the true 
confidence level of a confidence interval procedure is 80%. (respectively 
90%), then a 95% confidence interval for the fraction of the 300 replications 

that cover the true value is 0.80 ± 1.96 —1— (,2)(.8) = (.755, .845], 

oUU 

(respectively .9 ± 1.96 ^ ( . 1 ) ( . 9) = [.866, .934]). 

Table 2 shows the average length of confidence intervals for en P{D>t} 
for each t and each procedure. In parentheses under the average length is the 
estimated standard error of the average length. 

The conservative binomial procedure tends to overcover arid has the 
largest average confidence interval length. The maximum likelihood procedure 
has the correct coverage and smallest average length. The confidence 
intervals based on the asymptotic renewal estimate show the correct coverage 
for the times larger than t = 0.5 and 1.0. When the jackknife and bootstrap 
intervals do not cover they tend to be too high. The average length of the 
jackknife and bootstrap intervals is between those of the binomial and maximum 
likelihood intervals. The average length of the jackknife intervals is 
greater than those for the bootstrap, suggesting that the latter is preferred; 
however computational cost may be higher. 

Tables 3 and 4 present coverage results and average confidence interval 
lengths for a simulation experiment with N=50 individuals. The coverage 
results are similar to those of Table 1. As anticipated, the average 
confidence interval lengths tend to be smaller for N = 50 individuals than for 
N = 20 individuals. The average bootstrap interval length is closer to that of 
the MLE interval for times larger than 1 for N = 50 individuals than for N-20 
individuals . 

Tables 58 report results of a simulation experiment in which the sojourn 
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times in states 1 and 2 have the following distributions. The sojourn time in 



state 2 has an exponential distribution with mean — = 0 1 as before. The 



sojourn time in state 1 is the sum of two independent exponential random 
variables each having mean i. Thus, the mean sojourn time in each state is 
the same as for the model used in Tables 1 1. Values of P{D>t} for various 
values of t for this model are given in Appendix B. Estimators and confidence 
intervals for £n P{D> } are computed as before. In particular the maximum 
likelihood estimators assume that the sojourn times in state 1 have an 
exponential distribution with mean 1. The number of individuals is N = 20 for 
the simulation in Tables 5-6. There are N=50 individuals in the simulations 
of Tables 7-8. The coverage results contained in Tables 5 and 7 show that the 
maximum likelihood confidence intervals arc sensitive to underlying model 
assumptions. Incorrect model assumptions lead to both overcoverage and 
undercoverage. The confidence intervals based on the exponential 
approximation have the correct coverage for the times larger than 0.5. The 
average length of the jackknife and bootstrap intervals is smaller than those 
for the binomial. The average length of the bootstrap interval is smaller 
than those for the jackknife. 



4. CONCLUSIONS 

The simulation results indicate that the maximum likelihood estimator is 
most efficient when the correct model is used. As has been noted elsewhere, 
the MLE estimator is rather sensitive to incorrect model assumptions. The 
simulation results for the asymptotic exponential estimator are somewhat 
surprising in that it appears that the time t does not have to be very large 
in order for the asymptotic exponential estimator to work well. Both the 
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nonparametr ic bootstrap and jackknife confidence intervals for the asymptotic 
exponential estimator have the correct coverage for moderately large t. As 
implemented the bootstrap intervals tend to have a smaller average length than 
the jackknife. However, neither a complete jackknife nor a bootstrap with 
more than 100 replications has been done due to the massive amount of 
computation involved. As expected, the binomial estimator produces 
conservative confidence intervals which tend to overcover. However, the 
binomial confidence intervals are appropriate for all times and all 
semi-Markov models. 
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Appendix A 



Values of P{D>t} for the semi-Markov Model 
with exponential sojourn time in state 1 

0=0.5 



p i =1 



p 2 =10> 



time t 


P{D>t) 


0.5 


.7866 


1.0 


.6203 


1.5 


.4891 


2.0 


.3857 


3.0 


.2399 


4.0 


. 1492 
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Appendix B 



Values of P{D>t} for the semi-Markov Model 
with nonexponential sojourn time in state 1 

0=0.5 



Pj=2. 



p 2 =i°. 



p 3 -2. 



t ime t 


P(D>t> 


0.5 


.8652 


1.0 


.6743 


1.5 


.5157 


2.0 


.3930 


3.0 


.2278 


4.0 


. 1320 
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Table 1 



Coverage Results for Confidence Intervals for £n P{D>t} 
Exponential Sojourn Times 
N=20, 0=KD . 5 , p 1= 1.0, p 2 =10.0 



Time 


Conf idence 
Level 


Coverage 


BIN 


MLE 


JK 


BT 


0.5 


80% 


H 


14( . 05) 


24( .08) 


106( .35) 


140( .47) 






C 


272( .91 ) 


251( .84) 


174( .58) 


160( . 53) 






L 


14( .05) 


25( .08) 


20( .07) 


0(.0) 




90% 


H 


4(01) 


14( .05) 


78( .26) 


99( . 33) 






C 


282( .94) 


276 ( . 92) 


21 1 ( .70) 


201 (.67) 






L 


14( .05) 


10( .03) 


11( .04) 


0( .00) 


1.0 


80% 


H 


17( .06) 


24( .08) 


72( .24) 


66( . 22) 






C 


250( . 83) 


250( . 83) 


210( .70) 


224 ( .75) 






L 


33( . 11) 


26 ( . 09 ) 


18( .06) 


10( .03) 




90% 


H 


7 ( . 02 ) 


14( .05) 


46( . 15) 


33(. 11) 






C 


275 ( .92) 


277 ( .92) 


245( .82) 


254 ( .86) 






L 


18( .06) 


9( .03) 


9( . 03) 


8( .03 


1.5 


80% 


H 


15( .05) 


24( .08) 


48(. 16) 


40 ( . 13) 






C 


265( .88) 


250( .83) 


247( .82) 


236( .79) 






L 


20( .07) 


26 ( .09) 


5( .02) 


24( .08) 




90% 


H 


15( . 05) 


14 ( . 05) 


33(. 11) 


19( .06) 






C 


278 ( .93) 


277 ( .92) 


266( .89) 


266( .89) 






L 


7( .02) 


9( .03) 


1(.00) 


1 5 ( .05) 


2.0 


80% 


H 


1 1 ( . 04) 


24 ( . 08 ) 


50 ( . 17) 


27 ( . 09 ) 






C 


267 ( . 89 ) 


250( . 83) 


245 ( .83) 


240( . 80) 






L 


22( .07) 


26( .09) 


5( .02) 


33( . 11) 




90% 


H 


1 1 ( .04) 


14( .05) 


28 ( .09) 


10( . 03) 






C 


280( .93) 


277 ( .92) 


271 ( .90) 


267 ( .89) 






L 


9( .03) 


9( .03) 


1(.00) 


23( .08) 


3.0 


80% 


H 


23( .08) 


23( .08) 


40( . 13) 


19( .06) 






C 


27 1 ( .90) 


251 (.84) 


252( .84) 


238( .79) 






L 


6( .02) 


26( .09) 


8( .03) 


43(. 14) 




90% 


H 


8( .03) 


14( .05) 


22 ( .07) 


9( .03) 






C 


286( . 95) 


277 ( . 92) 


276 ( .92) 


260( . 87) 






L 


6( .02) 


9( . 03) 


2(.01) 


31 ( . 10) 


4.0 


80% 


H 


19( .06) 


23 ( . 08 ) 


36 ( . 12) 


15( .05) 






C 


27 1 ( . 90 ) 


251 ( .84) 


253( .84) 


231( .77) 






L 


10( .03) 


26( .09) 


11( .04) 


54( . 18) 




90% 


H 


4(.01) 


14( .05) 


21( .07) 


7( .03) 






C 


289 ( .96) 


277( .92) 


277( .92) 


261( .87) 






L 


7( .02) 


9( . 03) 


2( • 01 ) 


32( . 11) 
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Table 2 



Average Length of Confidence Intervals for £n P{D>t} 
Exponential Sojourn Times 
N=20, 0=0.5, p =1.0, p 2 =10.0 



Time Confidence BIN MLE JK BT 

Level 



0.5 


80% 


.370 


. 144 


.208 


. 185 




90% 


(.005) 

.460 


(.002) 
. 186 


( .005) 
.271 


( .003) 
.245 






(.007) 


( .002) 


( .006) 


(.005) 


1.0 


80% 


.545 


.288 


.307 


.290 




90% 


( .008) 
.680 


( .004) 
.370 


( .006) 
.400 


(.005) 

.382 






(.009) 


( .005) 


(.008) 


(.006) 


1 .5 


80% 


.716 


.431 


.467 


.436 




90% 


(.010) 

.893 


( .006) 
.553 


(.011) 

.608 


(.008) 

.579 






(.012) 


(.007) 


(.014) 


(.010) 


2.0 


80% 


.915 


.574 


.645 


.597 




90% 


(.016) 

1.14 


(.008) 

.737 


(.015) 

.840 


(Oil) 

.793 






(.020) 


(.010) 


(.019) 


(.015) 


3.0 


80% 


1.38 


.860 


1.02 


.932 




90% 


( .033) 
1.71 


(.011) 

1.10 


( .024) 
1.33 


(.018) 

1.24 






(.038) 


(0.15) 


( .032) 


(.025) 


4.0 


80% 


2.11 


1.15 


1.40 


1.27 




90% 


(.069) 

2.58 


(.015) 

1.47 


( .033) 
1.82 


( .026) 
1.70 






(.076) 


(.020) 


(.044) 


(.034) 
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Table 3 



Coverage Results for Confidence Intervals for £n P{D>t} 
Exponential Sojourn Times 
N=50, 0=0.5, Pj= 1.0, p 2 =10.0 



Time 


Conf idence 
Level 


Coverage BIN 


MLE 


JK 


BT 


0.5 


80% 


H 


17 ( .06) 


27 ( . 09 ) 


1 19( .40) 


15 1 ( .50) 






C 


260( .87) 


245( . 83) 


167 ( .56) 


149( . 50) 






L 


23( . 08) 


28( .09) 


1 4 ( .05) 


0(.00) 




90% 


H 


5 ( . 02 ) 


14( .05) 


83( .28) 


1 10( . 37) 






C 


288( .96) 


274( . 91 ) 


209( .70) 


190( . 63) 






L 


7( .02) 


12( .04) 


8( .03) 


0(.00) 


1.0 


80% 


H 


28 ( .09) 


27 ( . 09 ) 


70( .23) 


80( . 27) 






C 


245( .82) 


245 ( . 82 ) 


216( .72) 


213( .71) 






L 


27 ( .09) 


28( .09) 


14 ( . 05) 


7( .02) 




90% 


H 


8( .03) 


14( . 05) 


45(. 15) 


45(. 15) 






C 


282 ( .94) 


274( .91) 


274( . 82) 


250( .83) 






L 


10( .03) 


12( .04) 


8( . 03) 


5( .02) 


1.5 


80% 


H 


19( .06) 


27 ( . 09) 


53(. 17) 


46(. 15) 






C 


262 ( .87) 


245( .82) 


235( . 78) 


233 ( .78) 






L 


19( .06) 


28( .09) 


12( .04) 


2 1 ( . 07 ) 




90% 


H 


1 1 ( .04) 


1 4 ( .05) 


32( .11) 


23 ( .08) 






C 


278( .93) 


274( .91) 


262( .87) 


266( .89) 






L 


1 1 ( .04) 


12( .04) 


6( .02) 


11 (.04) 


2.0 


80% 


H 


18( .06) 


27 ( .09) 


47 ( . 16) 


31(. 10) 






C 


260( .87) 


245( .82) 


240( .80) 


243( .81) 






L 


22( .07) 


28( .09) 


13( .04) 


26 ( .09) 




90% 


H 


9( .03) 


14( .05) 


22 ( .07) 


15( .05) 






C 


277( .92) 


274( .91) 


270( .90) 


268( .89) 






L 


14( . 05) 


12( .04) 


8( .03) 


17 ( . 06) 


3.0 


80% 


H 


20 ( .07) 


27 ( .09) 


36 ( . 12) 


25( . 08) 






C 


262( .87) 


245( .82) 


250 ( .83) 


240( . 80) 






L 


18( .06) 


28 ( .09) 


14( . 05) 


35(. 12) 




90% 


H 


13( .04) 


14( .05) 


20( .07) 


11 (.04) 






C 


283 ( .94) 


275( .92) 


274 ( . 91 ) 


267(.89) 






L 


4(01) 


11( .04) 


6( .02) 


22( .07) 


4.0 


80% 


H 


14( .05) 


27( .09) 


35( . 12) 


20( . 07 ) 






C 


278( .93) 


245( . 82) 


248( .83) 


239 ( .80) 






L 


8( .03) 


28( .09) 


17 ( .06) 


41 (. 14) 




90% 


H 


6( . 02) 


14( . 05) 


19( .06) 


10( .03) 






C 


286 ( . 95 ) 


275( . 92) 


274( .91) 


262( . 87 ) 






L 


8( . 03) 


1 1 ( .04) 


7( .02) 


28( . 09) 
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Table 4 



Average Length of Confidence Intervals for £n P{D>t} 
Exponential Sojourn Times 
N=50. 0=0.5, p =1.0. P 2 =10.0 



Time 


Confidence BIN 
Level 


MLE 


JK 


BT 


0.5 


80% 


.216 


.089 


. 155 


. 127 






(.002) 


(.001) 


(.004) 


(.002) 




90% 


.271 


.114 


.205 


. 162 






(.002) 


(.001) 


(.006) 


(.003) 


1.0 


80% 


.319 


. 177 


.202 


. 181 






(.003) 


(.001) 


(.004) 


(.002) 




90% 


.401 


.227 


.268 


.234 






( .003) 


(.002) 


(.005) 


( .002) 


1.5 


80% 


.416 


.245 


.282 


.262 






(.003) 


(.002) 


(.005) 


( .003) 




90% 


.524 


.340 


.374 


.340 






(.004) 


(.003) 


(.007) 


(.004) 


2.0 


80% 


.524 


.353 


.377 


.352 






( .005) 


(.003) 


(.007) 


(.004) 




90% 


.628 


.453 


.500 


.458 






(.006) 


(.004) 


(.009) 


( .005) 


3.0 


80% 


.750 


.529 


.585 


.545 






( .008) 


(.004) 


(Oil) 


(.006) 




90% 


.940 


.579 


.775 


.706 






(.009) 


(.005) 


(.014) 


(.008) 


4.0 


80% 


1.05 


.705 


.799 


.738 






(.015) 


(.006) 


(.015) 


(.009) 




90% 


1.31 


.905 


1.06 


.959 






(.019) 


(.007) 


(.020) 


(.011) 
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Table 5 



Coverage Results for 'Confidence Intervals for £n P{D>t} 
Nonexponential Sojourn Times 
N=20, 0=0.5, Pj=2. 0, p =10.0. p 3 =2.0 



Time 


Conf idence 
Level 


Coverage 


BIN 


MLE 

(Exp) 


JK 


BT 


0.5 


80% 


H 


12( . 04) 


0(.0) 


105( .35) 


135( .45) 






C 


278( . 93) 


53 ( . 18) 


183( .61) 


162( . 54) 






L 


10( . 03) 


247( .82) 


12( . 04) 


3( .00) 




90% 


H 


0( .00) 


0( .00) 


73( .24) 


82( .37) 






C 


290( .97) 


109( .36) 


220 ( .73) 


217 ( . 72) 






L 


10( .03) 


191 ( .64) 


7( .02) 


1(0) 


1.0% 


80% 


H 


16( .05) 


1(0) 


45 ( . 15) 


32( .11) 






C 


254 ( .85) 


221 ( .74) 


240 ( .80) 


241 (.80) 






L 


30 ( . 10) 


78( .26) 


1 5 ( .05) 


27 ( .09) 




90% 


H 


5( .02) 


0(.0) 


29 ( . 10) 


11( .04) 






C 


285( .95) 


266 ( . 89 ) 


263( .88) 


271 ( .90) 






L 


10( .03) 


34( .11) 


8( .02) 


18( .06) 


1.5 


80% 


H 


21 ( .07) 


8( .03) 


40 ( . 13) 


23( .08) 






C 


260( .87) 


257( .86) 


244( .81) 


232( .77) 






L 


19( .06) 


35(. 12) 


16( .05) 


45 ( . 15) 




90% 


H 


5( .02) 


2 ( .01) 


24( .08) 


7( .02) 






C 


276( .92) 


283( .94) 


273 ( .91) 


264( .88) 






L 


19( .06) 


1 5 ( . 05 ) 


3(01) 


29 ( . 10) 


2.0 


80% 


H 


1 1 ( .04) 


17( . 06) 


39 ( . 13) 


18( .06) 






C 


270( .90) 


260 ( .87) 


246( .82) 


232( .77) 






L 


19( .06) 


23 ( . 08 ) 


15( . 05) 


50( . 17) 




90% 


H 


11( .04) 


7( .02) 


19( .06) 


7( .02) 






C 


284( .95) 


284( .95) 


277( .92) 


263( .88) 






L 


5( .02) 


9( .03) 


4( • 01 ) 


30(. 10) 


3.0 


80% 


H 


13 ( .04) 


31 ( . 10) 


37 ( . 12) 


17( .06) 






C 


273 ( .91) 


254( .85) 


251( .84) 


230( .77) 






L 


14( .05) 


15( .05) 


12( .04) 


53 ( . 18) 




90% 


H 


5( .02) 


17( .06) 


19( .06) 


6( .02) 






C 


281 ( .94) 


279( .93) 


277( .92) 


261( .87) 






L 


14( .05) 


4(01) 


4(01) 


33( . 11) 


4.0 


80% 


H 


4( .01 ) 


44 ( . 15) 


38( . 13) 


17( . 06) 






C 


279 ( .93) 


245( .82) 


251 ( . 84) 


229( .76) 






L 


17( . 06) 


11( .03) 


1 1 ( .04) 


54 ( . 18) 




90% 


H 


4(.01) 


23( .08) 


19( .06) 


6( .02) 






C 


296( .98) 


273( .91) 


277 ( . 92 ) 


260 ( .87) 






L 


0( .00) 


4( 01 ) 


4(01) 


34 ( -11) 
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Table 6 



Average Length of Confidence Intervals for £n P{D>t} 
Nonexponential Sojourn Times 
N=20. 0=0.5, p 1 =2 . 0 , p 2 =10.0 p 3 =2.0 



Time 


Confidence BIN 


MLE 


JK 


BT 




Level 




(Exp) 






0.5 


80% 


.289 


. 144 


. 141 


. 130 






(.004) 


(.002) 


(.003) 


(.003) 




90% 


.359 


.185 


. 183 


.175 






(.005) 


(.002) 


(.004) 


(.004) 


1.0 


80% 


.485 


.286 


.269 


.253 






(.007) 


( .0Q3) 


(.006) 


(.004) 




90% 


.604 


.367 


.351 


.331 






(.008) 


(.004) 


(.008) 


(.006) 


1.5 


80% 


.681 


.428 


.451 


.419 






(.010) 


(.005) 


(.010) 


(.007) 




90% 


.850 


.550 


.588 


.552 






(.013) 


(.007) 


(.014) 


(.010) 


2.0 


80% 


.884 


.571 


.643 


.594 






(.013) 


(.007) 


(.015) 


(.017) 




90% 


1.10 


.733 


.838 


.786 






(.017) 


(.009) 


(.019) 


(.015) 


3.0 


80% 


1.46 


.855 


1.03 


.953 






(.034) 


(.010) 


(.024) 


(.018) 




90% 


1.81 


1.10 


1.35 


1.26 






(.043) 


(.013) 


(.032) 


(.025) 


4.0 


80% 


2.34 


1.14 


1.43 


1.31 






(.079) 


(.013) 


(.034) 


(.025) 




90% 


2.84 


1.46 


1.86 


1.74 






(.084) 


(.017) 


(.044) 


(.035) 
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Table 7 



Coverage Results for Confidence Intervals for Cn P{D>t} 
Nonexponential Sojourn Times 
N=50, 0=0.5. p 1 =2.0, p 2 =10.0. p 3 = 2.0 



Time 


Conf idence 
Level 


Coverage 


BIN 


MLE 

(Exp) 


JK 


BT 


0.5 


80% 


H 


18( . 06) 


0(.0) 


130( .43) 


149( .50) 






C 


261 ( .87) 


3( • 01 ) 


165( .55) 


150( .50) 






L 


21( .07) 


297 ( .99) 


5( .02) 


1(0) 




90% 


H 


10( .03) 


0( .00) 


99( . 33) 


106( .35) 






C 


283( .94) 


6( .02) 


200( .67) 


194( .65) 






L 


7( .02) 


294( .98) 


1(.00) 


0(.0) 


1.0 


80% 


H 


22( .07) 


2(01) 


42( . 14) 


39( . 13) 






C 


260( . 87) 


145( . 48) 


242( .81) 


240( . 80) 






L 


18( . 06) 


153( .51 ) 


16( . 05) 


21 ( .07) 




90% 


H 


1 1 ( .04) 


0( .00) 


28( .09) 


21 ( .07) 






C 


280 (.93) 


213( .71 ) 


263 ( .88) 


265( . 88) 






L 


9( .03) 


87( . 29) 


9( .03) 


14( .05) 


1.5 


80% 


H 


27 ( .09) 


6( .02) 


31 ( . 10) 


25( .08) 






C 


254( .85) 


242( .81) 


247 ( . 82 ) 


243( .81) 






L 


19( . 06) 


52 ( . 17) 


22( .07) 


32( .11) 




90% 


H 


1 1 ( .04) 


4( • 01 ) 


21 (.07) 


19( .06) 






C 


276( . 92) 


276( .93) 


27 1 ( . 90 ) 


262( .87) 






L 


13( .04) 


20 ( .07) 


8( .02) 


19( .06) 


2.0 


80% 


H 


25( . 08) 


13( .04) 


29 ( . 10) 


24( .08) 






C 


253( . 84) 


264( .88) 


250( . 83) 


239( . 80) 






L 


22( .07) 


23( .08) 


2 1 ( . 07 ) 


37 ( . 12) 




90% 


H 


16( . 05) 


7( .02) 


20( . 07 ) 


12( .04) 






C 


274( .91) 


282( .94) 


269( .90) 


264( .88) 






L 


10( .03) 


1 1 ( .04) 


1 1 ( .04) 


24( .08) 


3.0 


80% 


H 


25( . 08) 


34 ( 11) 


31 ( . 10) 


21 ( .07) 






C 


247 ( . 82 ) 


253( . 84) 


248( .83) 


241 ( .80) 






L 


28( . 09) 


1 3 ( .04) 


21 { .07) 


38( . 13) 




90% 


H 


17( .06) 


20( .07) 


1 8 ( .06) 


1 1 ( .04) 






C 


273( .91) 


275 ( .92) 


271( .90) 


266( . 89) 






L 


10( . 03) 


5( .02) 


1 1 { .04) 


23( .08) 


4.0 


80% 


H 


1 1 ( .04) 


47(. 16) 


31 ( . 10) 


20( . 07 ) 






C 


266( .89) 


246( . 82) 


244( .81) 


238( . 79) 






L 


23( . 08) 


7(.02) 


25( .08) 


42 ( . 14) 




90% 


H 


6( . 02) 


29 ( . 10) 


19( .06) 


10( .03) 






C 


285( .95) 


267( .89) 


270( .90) 


264( .88) 






L 


9( .03) 


4(01) 


11( .04) 


26( .09) 
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Table 8 



Average Length of Confidence Intervals for in P{D>t} 
Nonexponential Sojourn Times 
N=50, 0=0.5, p^.O, p =10.0, p 3 =2.0 



Time Confidence BIN MLE JK BT 

Level 



0.5 


80 % 


. 169 


.089 


.091 


.081 




90 % 


(. 002 ) 

.211 


(. 001 ) 

.114 


(. 002 ) 
. 121 


(. 001 ) 

.104 






(. 002 ) 


(. 001 ) 


(. 003 ) 


( . 002 ) 


1.0 


80 % 


.284 


. 176 


. 167 


. 153 




90 % 


(. 002 ) 

.357 


(. 001 ) 

.227 


( . 003 ) 
.221 


(. 002 ) 

.199 






(. 003 ) 


( . 002 ) 


(. 004 ) 


( . 002 ) 


1.5 


80 % 


.356 


.264 


.269 


.250 




90 % 


(. 003 ) 

.498 


(. 002 ) 

.339 


(. 005 ) 

.357 


( . 003 ) 
.324 






(. 004 ) 


(. 002 ) 


(. 006 ) 


(. 004 ) 


2.0 


80 % 


.514 


.352 


.378 


.351 




90 % 


(. 005 ) 

.646 


(. 003 ) 

.452 


(. 007 ) 

.501 


(. 004 ) 

.456 






(. 006 ) 


(. 003 ) 


(. 010 ) 


(. 005 ) 


3.0 


80 % 


.784 


.527 


.600 


.558 




90 % 


(. 009 ) 

.983 


(. 004 ) 

.677 


( Oil ) 

.795 


(. 006 ) 

.726 






( Oil ) 


(. 005 ) 


(. 014 ) 


(. 009 ) 


4.0 


80 % 


1.16 


.703 


.824 


.765 




90 % 


(. 019 ) 

1.45 


(. 063 ) 

.903 


(. 015 ) 

1.09 


(. 009 ) 

.998 






(. 024 ) 


(. 064 ) 


(. 020 ) 


(. 012 ) 
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